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Abstract. We consider the production of 2.22 MeV ra- 
diation resulting from the capture of neutrons in the at- 
mosphere of the secondary in an X-ray binary system, 
where the neutrons are produced in the accretion disk 
around the compact primary star and radiated in all di- 
rections. We have considered several accretion disk models 
(ADAF, ADIOS, SLE, Uniform-Temperature) and a var- 
ity of parameters (accretion rate, mass of the compact ob- 
ject, mass, temperature and composition of the secondary 
star, distance between the two objects, etc.). The neu- 
tron rates are calculated by a network of nuclear reactions 
in the accretion disk, and this is handled by a reaction- 
rate formulation taking into account the structure equa- 
tions given by each accretion model. The processes under- 
gone by the neutrons in the atmosphere of the companion 
star are studied in great detail, including thermalization, 
elastic and inelastic scatterings, absorption, escape from 
the surface, decay, and capture by protons. The radiative 
transfer of the 2.22 MeV photons is treated separately, 
taking into consideration the composition and density of 
the star's atmosphere. 

The final flux of the 2.22 MeV radiation that can be 
detected from earth is calculated taking into account the 
distance to the source, the direction of observation with 
respect to the binary system frame, and the rotation of the 
source, as this can lead to an observable periodicity in the 
flux. We produce phasograms of the 2.22 MeV intensity 
as well as spectra of the line, where rotational Doppler 
shift effects can lead to changes in the spectra that are 
measurable by INTEGRAL'S spectrometer (SPI). 

Key words: X-rays: binaries - Accretion, accretion disks 
- Gamma-ray: theory - Lines: profiles - Nuclear reactions, 
nucleosynthesis, abundances 



1. Introduction 

The capture of a neutron by a proton and the subsequent 
emission of a photon at 2.223 MeV has been the object of 



much interest since the dawn of gamma-ray astrophysics 
(Fichtel & Trombka 1981) because it represents a poten- 
tially important spectroscopic tool in the investigation of 
high-energy environments such as solar flares, neutron star 
surfaces or accretion disks around compact objects. Sev- 
eral studies have been conducted both from theory and 
observation, realizing that this line is the most promising 
candidate of all nuclear gamma-ray processes. 

Solar flares have had the lion's share of works because 
of the obvious possibilities of detailed investigation: Hua 
& Lingenfelter (1987); Murphy et al. (1991); Share & Mur- 
phy (1995); Ramaty et al. (1995); Ramaty, Mandzhavidze, 
Kozlovzky (1996); and others. Other sources have also 
been investigated, although no detection has been con- 
firmed and established to date, despite several previous 
announcements (Jacobson 1982; McConnell et al. 1997). 

On the theoretical front, Guessoum & Dermer (1988) 
conducted a detailed theoretical investigation in the con- 
text of compact binary sources, taking Cygnus X-l as a 
case-study and pointing out the possibility of emission of 
a very narrow 2.22 MeV line from the atmosphere of the 
companion, while Aharonian & Sunyaev (1984) had con- 
sidered a two-temperature accretion disk and considered 
the possibility of emission from the disk itself, where the 
line would be very broad. Bildsten and co-workers consid- 
ered the possibility of emission at the surface of a neu- 
tron star, as a result of bombardment by accreted pro- 
tons (Bildsten 1991; Bildsten, Salpeter, and Wasserman 
1993). In the same framework, Bykov et al. (1999) recently 
considered the possibility of detecting nuclear gamma-ray 
lines from accreting objects with INTEGRAL. Vestrand 
(1989) investigated the 2.22 MeV emission resulting from 
very high energy (E > 10 12 eV) protons bombarding com- 
panion stars of sources such as Cyg X-3, Vel X-l, and 
Her X-l, and predicted significant fluxes (« 10~ 4 pho- 
tons cm -2 s _1 ). Finally, Guessoum & Kazanas (1999), in- 
terested in Lithium production by neutron bombardment 
of X-ray binary companions' atmospheres and using the 
ADAF accretion disk model, estimated the flux of 2.22 
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MeV line emission from a source at about 1 kpc and found 
it generally very low ( 10~ 7 photons cm~ 2 s _1 ). 

On the observational front, extensive searches of the 
line have been conducted with the data from the So- 
lar Maximum Mission (Harris & Share 1991) and the 
COMPTEL instrument of the Compton-GRO mission 
(McConnell et al. 1997). The latter showed an excess emis- 
sion at (l,b) = 300°,-30° at about 3.7 a. No obvious coun- 
terparts (X-ray binaries or the like) could be linked to that 
general direction, although a catalysmic variable has been 
suggested as a likely candidate for that emission. Van Dijk 
(1996) has compiled a list of upper-limit 2.22 MeV binary- 
source radiation fluxes based on COMPTEL observations 
of 27 black-hole candidates. These upper limits range from 
1 to 5 xlCT 5 photons cm~ 2 s _1 . 

The production of neutrons in the accretion disks 
around the primaries of X-ray binaries is a direct result of 
the high temperatures that prevail in such environments. 
The gravitational attraction, combined with the viscous 
dissipation, leads to a significant heating of the plasma, 
particularly its nuclei, resulting in temperatures as high 
as 10 10 - 10 12 K in the inner regions (r « 1 — 100 R5, 
where Rs = 2GM/c 2 is the Schwarzchild radius of the 
central compact object). Simple considerations show that 
densities of the accretion plasma can range between 10 12 
and 10 18 cm 3 depending on the geometry and the model, 
which implies substantial nuclear reactions. Since the ac- 
creted material is expected to contained He and/or metals 
(depending on whether the secondary star's composition is 
normal or of Wolf-Rayet type), breakup reactions will pro- 
duce significant amounts of neutrons. Part of the present 
paper is devoted to computing precisely such amounts un- 
der various assumptions: different accretion disk models, 
different initial compositions of the accreted material, dif- 
ferent temperature profiles, etc. This will be presented in 
detail in Section 2. 

There are several reasons underlying our interest in 
2.22 MeV emission from the secondary stars of X-ray bi- 
naries. First and foremost, as mentioned above the CGRO 
2.22 MeV map showed that there are real possibilities 
for the detection of this line, which would consitute the 
first observation of any nuclear gamma-ray line outside of 
the solar system, and this prospect is heightened by the 
upcoming INTEGRAL mission's order-of-magnitude im- 
provement in both flux sensitivity and energy resolution, 
which become crucial when most predictions of line fluxes 
from such sources are at the 10~ 5 photons cm~ 2 s _1 level, 
that is below the sensitivity of CGRO and right at the de- 
tection threshold of INTEGRAL. Secondly, the accretion 
disk problem remains fully unresolved (very little progress 
has been made on the determination of the disk structure, 
much less the identification of the viscosity mechanism), 
and since our results will show significant differences in 
fluxes and therefore in the possibility for detection of this 
line emission, this may constitute further testing of mod- 
els. Furthermore, we will find that line fluxes depend sub- 



stantially on the secondary star's characteristics (elemen- 
tal composition in particular), the eventual detection of 
such a signal would then also constitute a diagnostic of 
the X-ray binary's secondary (e.g. convection and mixing 
of the gas). 

This paper is structured as follows: in Sect. 2, the treat- 
ment of neutron production in the accretion disk (and 
its subsequent escape) is presented for various models; in 
Sect. 3, the detailed treatment of neutron propagation, 
slow-down, and capture in the atmosphere of the sec- 
ondary is shown; in Sect. 4 we present the results of line 
fluxes when everything is combined: neutron production 
and escape from the disk, various interactions and effects 
(geometrical as well as physical) in the atmosphere of the 
secondary, and finally capture and gamma-ray emission, 
with star rotation, photon energy shift, etc.; in the final 
section we discuss our results with INTEGRAL observa- 
tion prospects in mind and point to future work on the 
subject. 

2. Neutron Production in the Accretion Disk 

Neutrons are produced by means of nuclear breakup re- 
actions involving Hydrogen and Helium nuclei. Once pro- 
duced, the neutrons can escape or carry part of the angular 
momentum outward by collisions with the inf ailing nuclei, 
thereby participating in the viscous dissipation, as they 
are not affected by the presence of any magnetic fields. 
In fact, the escape of the neutrons from the gravitational 
attraction of the compact object can be computed quite 
accurately in the assumption of thermal conditions (Aha- 
ronian & Sunyaev 1984; Guessoum & Kazanas 1990). So, 
at least in principle, the calculation of the flux of neu- 
trons irradiating the surface of the companion star is quite 
feasible. In practice, however, there are a number of fac- 
tors which intervene and make the task somewhat more 
complicated: the temperature and density distribution of 
matter in the accretion disk is model-dependent, the com- 
position and nuclear abundances of the plasma depends 
on stellar evolution conditions, and cross sections for the 
various reactions are not always known accurately (and 
this depends also on the extent of the nuclear reactions 
network taken into account in the computation). 

Since the nuclear reactions relevant to our problem 
and to the production of gamma rays require energies of 
about 10 MeV per nucleus, our interest is limited to the 
inner regions of the accretion flow, i.e. regions around the 
compact object between about 100 Schwarzschild radii) 
and either the surface of the neutron star or the horizon 
of the black hole, depending on the nature of the compact 
object. For general purposes, we assume a mass of 1 Mq 
for the compact object, unless a specific case is considered. 

In this work we consider a number of relatively simple 
accretion disk models: 

— 1 - The Advection-Dominated Accretion Flow (ADAF) 
model proposed by Narayan and co-workers (Narayan 
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& Yi 1994, Chen et al. 1995, Narayan, Yi, and Mahade- 
van 1995, Narayan & Yi 1995), which has the added 
appeal that its temperature and density function arc 
very simple and analytical expressions (Yi & Narayan 
1997): 

Tj = 1.11 x 10 12 r _1 K = 95.68 r _1 MeV (1) 

where r is the dimensionless radial variable r = R/Rs; 

n = 6.4 x 10 18 a" 1 to" 1 mr~ 3/2 cuT 3 (2) 

where to is the dimensionless mass measured in units 
of a solar mass (m = M/M ), m is the dimensionless 
accretion rate measured in units of the Eddington ac- 
cretion rate m = M/MEdd = M/1.4 x 10 17 m g/s (in 
this work we will consider values of M in the range 
of 1(T 10 M Q /yr - 1CT 8 M Q /yr), taking (3 (the ratio 
of gas pressure to total pressure) to equal 1/2 (as in 
Yi and Narayan 1997); the central factor a is taken to 
equal either 0.3 or 0.1 (following Narayan & Yi 1994). 
Furthermore, it must be noted that in this model the 
value of the electron temperature does not affect our 
results, which depend mainly on the profile of the ion 
temperature T^(r). 

2 - The Advection-Dominated Inflow-Outflow Solu- 
tions proposed by Blandford and Begelman (1998), 
which consists of a "generalized" ADAF model where 
only a small fraction of the gas initially supplied ends 
up falling onto the compact object, and where the en- 
ergy released is transported radially outward and ef- 
fectively becomes a wind, driving away the rest of the 
accreted material. 

In this family of solutions, a central assumption 
is the radial dependence of the accretion rate as 
m(r) oc r p , < p < 1 . In the present work, we take the 
"gasdynamical wind" (case iv of Blandford & Begel- 
man 1998) solution, i.e. p = 1/2, e = 1/2, which leads 
to a density function n oc r ; more precisely: 



n : 



„ ,„ ,„iq i (4-52 x 10 8 ) x M , 
2.18 x 10 19 a- 1 ± '- cm" 3 



mr 



(3) 



Note that this special form of the density function also 
corresponds to the solution obtained by Kazanas, Hua, 
& Titarchuk (1997) in fitting the time lags in similar, 
compact X-ray accreting sources. 

Another important equation is that of the disk height 
h (in units of Rg), which in the particular solution 
adopted here, is given by: 



h = 0.447 x 3mr; 

Also, the radial velocity is 

v r = 0.52 a r~2 c cm/ s . 



(4) 



(5) 



Finally, the ion temperature profile is the same as in 
ADAF. 



3 - The Two- Temperature Shapiro-Lightman-Eardley 
(1976) Model, for which the expressions are given by: 

T % = 5 x 10 11 M~* M? q-5 4>i r~ f K (6) 

T e = 7x 10 8 M? M,"'a-5f^,' K (7) 

h = 10 5 MP M^* a-T2 (j)T2 rl cm (8) 

p = 5 x 10~ 5 M~* ai <j>~* r* S g/cm 3 , (9) 

where = M/3; = l-y/3R s /R= l-y/3/r;M* = 
M/10 17 g/s;n = 2R/R S = 2r . 

4 - The uniform-temperature model, similar to the 
one used in the Guessoum & Kazanas (1990) neutron- 
viscosity model, where Tj is set by a viscosity condition 
but is taken to be uniform in the disk, and T e is chosen 
arbitrarily (usually between 0.1 and 0.5 MeV); the den- 
sity is then obtained from an energy balance require- 
ment (viscous heating of the ions equal to Coulomb 
energy transfer from the ions to the electrons). In the 
Guessoum & Kazanas model, the following relations 
were derived: 



3.375 x 10 25 to 



m 3 fth 



(10) 



where fth is the ion-electron Coulomb energy transfer 

function (Dermer 1986), 

and 



47r 3 n 

tvisc — ^ 

3 TO 

where mj is the mass of the ion. 



(11) 



Finally we note that in each case the system is taken to 
be in steady-state, whereby the accretion and in-fall are 
constant in time, and thus the neutron production and 
gamma-ray emission are steady. 

2.1. The Nuclear Reactions 

Assuming an initial composition for the accreting plasma, 
that is some initial fractions of Hydrogen and Helium (90% 
H, 10% He in the "normal" case, and 10% H, 90% He in 
the "helium-rich case" , which would correspond to a Wolf- 
Rayet secondary star companion), we compute the rates 
of the main nuclear reactions which the ions (the protons, 
alphas, 3 He and 2 H nuclei) can undergo: 



p + a 



i He + 2 H 



a + a ► 7 Li + p 

p + 3 He — ► 3p + n 
a + a — ► 7 Be + n 
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a + 3 He > a + 2p + n 

p + 2 H — >2p + n 

The first two reactions contributing to the destruc- 
tion of Helium (our source of neutrons), and the last 
four reactions are the actual neutron-producing processes. 
We disregard all high-energy neutron-production from 
proton-proton collisions. The cross sections are essen- 
tially the same as those used in Guessoum & Kazanas 
(1999), and the reaction rates are calculated through 
the usual non-relativistic expression for binary processes 
(jij = niTij/{\ + Sij) < (JijVij >, the averaging here is 
performed over thermal distributions for the ions). 

We then let the plasma evolve over the dynamical 
timescale td, with its temperature and density set by the 
disk structure equations given above for each model, as 
the material sinks in the gravitational well. We use a sim- 
ple numerical scheme of explicit finite-differencing to fol- 
low the abundances of the various species (the neutrons 
in particular). This calculation is performed for various 
values of the model parameters: M — 10 -10 , 10~ 9 and 
1CT 8 Af Q /yr; a = 0.3 and 0.1; M = 1; etc. Note that since 
we are only interested in neutron production in the disk 
and not in any radiative processes, the electron temper- 
ature T e , whenever needed, is taken as a free parameter, 
usually equal to either 100 or 500 keV. 

The escape of the neutrons is handled in the same 
way as in Aharonian & Sunyaev (1984) or Guessoum 
& Kazanas (1990), that is by computing the fraction 
of neutrons that have kinetic energies sufficient to over- 
come the gravitational binding to the compact object: 
l/2m n (vfi ow + vthermai) 2 > GMm n /R; this translates 
into a fraction of escaping neutrons given by Eqs. (13) - 
(15) of Guessoum & Kazanas (1990). 

One final important aspect of the neutrons produced 
is their energy distribution. This aspect was addressed 
briefly in Guessoum & Kazanas (1999), where it was con- 
cluded that a thermal distribution of the neutrons with a 
temperature equal to that of the ions in the disk is reason- 
able. This point proves to be important not only for the 
escape fraction, but also for the interaction of those neu- 
trons that reach the secondary star's atmosphere, since 
their various interactions there depend strongly on their 
kinetic energies. 

2.2. Neutron Production Results in Various Disks 

Before presenting results of neutron production for the 
various disk models, we first show the main characteris- 
tics of these models. Figures [l] show the plasma density n 
(in units of 10 15 cm' 3 ), the ion temperature Tj (in units of 
MeV), the escape fraction of neutrons f esc , and the abun- 
dance of neutrons achieved in the disk, all as a function of 
r = R / R s , in the ADAF, ADIOS, SLE, and Uniform-T 4 
(10 and 30 MeV), respectively, for disks with solar com- 
position (90% H and 10% He) and an accretion rate of 



Table 1. Neutron Production Rate (in neutrons/s) in the 
ADAF Model, for initial (a) 90% H and 10% He and (b) 
10% H and 90% He compositions. 



M(M Q /yr) 






a 












b 








a 




0.3 


a 




0.1 


a 




0.3 


a 




0.1 


1Q -io 


1.1 


X 


10 37 


5.6 


X 


10 37 


6.6 


X 


10 37 


2.7 


X 


10 38 


10" 9 


4.4 


X 


10 38 


1.4 


X 


10 39 


2.6 


X 


10 39 


6.0 


X 


10 39 


icn 8 


9.1 


X 


10 39 


1.6 


X 


10 40 


3.9 


X 


10 40 


7.3 


X 


10 40 



Table 2. Neutron Production Rate in the ADIOS Model, 
for initial (a) 90% H and 10% He and (b) 10% H and 90% 
He compositions. 



M(M Q /yr) 






a 












b 








a 




0.3 


a 




0.1 


a 




0.3 


a 




0.1 


1Q -io 


1.7 


X 


10 37 


1.4 


X 


10 38 


1.4 


X 


10 38 


7.9 


X 


10 38 


10" 9 


1.6 


X 


10 39 


7.8 


X 


10 39 


9.9 


X 


10 39 


3.5 


X 


10 40 


10" 8 


8.3 


X 


10 40 


2.3 


X 


10 41 


3.6 


X 


10 41 


9.8 


X 


10 41 



Table 3. Neutron Production Rate in the SLE Model, for 
initial (a) 90% H and 10% He and (b) 10% H and 90% He 
compositions. 



M(M Q /yr) 








a 












b 












0.3 


a 




0.1 


a 




0.3 


a 




0.1 


10 -io 


3.6 


X 


10 23 


5.0 


X 


10 33 


6.5 


X 


10 24 


8.6 


X 


10 34 


10" 9 


4.2 


X 


10 36 


2.9 


X 


10 38 


5.0 


X 


10 37 


1.8 


X 


10 39 


10" 8 


5.0 


X 


10 39 


1.9 


X 


10 40 


3.3 


X 


10 40 


1.3 


X 


10 41 



10" 9 M /yr; for ADAF, ADIOS, and SLE the value of a 
in the figures is 0.1, and for the uniform-T^ models the 
electron-temperature used for the figures is 0.5 MeV. 

Table 1 presents the results of these calculations, show- 
ing the rates of neutron production in the ADAF-model 
case for various values of the viscosity parameter a and 
the mass accretion rate M (la corresponds to a "nor- 
mal" composition; lb corresponds to the "Helium-rich" 
case) . All the results shown correspond to M = 1 M Q and 
R w 100 Rs ~ 10 7 cm, though the actual results are 
mostly insensitive to the value of the outer radius of the 
disk, as most of the nuclear reactions and neutron produc- 
tion take place in the inner part (less than about 10 Rs) 
of the accretion flow where the ion temperature is several 
tens of MeV. 

Table 2 is equivalent to Table 1 in the ADIOS model. 
Table 3 corresponds to the SLE model; Tables 4 and 5 
correspond to the uniform-Ti model, with T^ = 30 MeV 
and T l = 10 MeV. Respectively. 

It is interesting to note that the resulting neutron 
fluxes depend strongly on the models considered, but also 
are quite sensitive to the 3 disk parameters considered, 
namely a, M, and the initial composition of the plasma. 



P. Jean & N. Guessoum: Neutron-Capture and 2.22 MeV emission 



10" 
10 3 
io 2 

10 1 
10" 

lir 1 

10" 2 

io- 3 
io 4 

10 4 
10 3 
IO 2 
IO 1 
10" 
IO" 1 
L0' 2 
IO" 3 

nr 4 



-T (MeV) 

n(10 15 cm 3 ) 
f 



20 



60 



SLE 



-T. (MeV) 

n(10 I5 cm 3 ) 
f 



10 4 
IO 5 
IO 2 
IO 1 
10" 
IO" 1 
IO" 2 
IO" 3 

io 4 

10 4 
10 3 
IO 2 
IO 1 
10" 
IO" 1 
IO" 2 

io 3 
ur 4 



20 



-T. (MeV) 

n(10 15 cm" 3 ) 

-f 



60 



\ 

\ 

r \ 


Uniform-T.(lOMeV) 






-> 




n (IO 15 cm -3 ) 

f ic 

a 



IO 4 
IO 3 
IO 2 
IO 1 
10" 
IO 1 
IO" 2 
IO" 3 

ur 4 



\ 

\ 


Uniform-T. (30 MeV) 






/ 


T (MeV) 

n(10 15 cm- 3 ) 




f 




a 





20 



Fig. 1. Plasma characteristics (density, ion temperature, escape fraction of neutrons, and abundance of neutrons) as 
a function of r, in the ADAF, ADIOS, SLE, and Uniform-T, (10 and 30 MeV), assuming M = 10~ 9 M Q /yr, a = 0.1, 
and T e = 0.1 MeV. 



3. Propagation, slow-down and capture of 
neutrons in the atmosphere of the secondary 

A fraction of the neutron flux emitted by the accretion 
disk irradiates the atmosphere of the secondary star. Some 
of these neutrons get thermalized; some are captured by 
nuclei, others decay or escape the secondary if after several 
scatterings their kinetic energies become larger than the 



gravitational potential energy binding them to the star. 
The 2.22 MeV photons resulting from the capture of neu- 
trons by protons escape the secondary atmosphere if they 
are not absorbed or Compton scattered in the surrounding 
gas. Consequently, the probability of escape of 2.22 MeV 
photons depends on the depth of their creation site (in 
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Table 4. Neutron Production Rate in the Uniform Ion 
Temperature Model with kT* = 30 MeV , for initial (a) 
90% H and 10% He and (b) 10% H and 90% He composi- 
tions. 



M(M Q /yr) 


(1 




b 




kT e = 


kT e = 


kT e = 


kT e = 




0.5MeV 


O.lMeV 


0.5MeV 


O.lMeV 


1(T 10 


3.2 x 10 38 


4.8 x 10 37 


7.4 x 10 38 


1.0 x 10 38 


1Q -9 


2.4 x 10 39 


4.3 x 10 38 


6.0 x 10 39 


1.0 x 10 39 


1Q -8 


1.6 x 10 40 


5.1 x 10 39 


4.4 x 10 40 


3.1 x 10 40 



Table 5. Neutron Production Rate in the Uniform Ion 
Temperature Model with kT = 10 MeV , for initial (a) 
90% H and 10% He and (b) 10% H and 90% He composi- 
tions. 



M(M e /yr) 




a 




b 




kT e = 


kT e = 


kT e = 


kT e = 




0.5MeV 


O.lMeV 


0.5MeV 


O.lMeV 


1Q -io 


2.0 x 10 38 


2.1 x 10 37 


4.0 x 10 38 


4.1 x 10 37 


10 -9 


1.6 x 10 39 


2.1 x 10 38 


3.6 x 10 39 


4.8 x 10 38 


10" 8 


1.1 x 10 40 


1.8 x 10 39 


2.7 x 10 40 


8.0 x 10 39 



the atmosphere) and on their direction of emission with 
respect to the surface. 

These processes have all been modeled, in the aim of 
determining the rate of 2.22 MeV radiation emitted by the 
secondary. The transport of neutrons in the secondary has 
been simulated using the code GEANT/GCALOR, which 
takes into account elastic and inelastic interactions as well 
as the neutrons' decay and capture by nuclei. The radia- 
tive transfer of 2.22 MeV photons in the atmosphere and 
the fraction of neutrons falling back onto the secondary 
are then modeled separately, using the results of the pre- 
vious simulations. The total emissivity at 2.22 MeV is es- 
timated by integrating over the whole secondary surface 
irradiated by neutrons. In the following paragraphs details 
of the method are presented. 

The secondary star has a mass and a radius denoted 
hereafter by M s and R s . (There should be no confusion 
between the Schwarzchild radius (Rs) of the primary and 
the radius (R s ) of the secondary star.) The star's atmo- 
sphere is characterized by a temperature T s , a density p s 
and an abundance in H and in He of [H] and [He] , respec- 
tively. The isotopic abundance of these elements are taken 
from Anders and Grevesse (1989). We have disregarded all 
other elements in the atmosphere. The distance between 
the accretion disk and the center of the secondary (also 
called separation) is parametrized by d s . This distance d s 
and the radius of the star determine the fraction of the 
neutron flux emitted by the accretion disk that irradiates 
the secondary. The secondary's mass and radius allow us 
to calculate the gravitational potential and therefore the 
fraction of escaping neutrons, which also depends on their 



energy and direction. The temperature, the density, and 
the composition of the atmosphere are parameters for the 
simulation of neutron tracks in the secondary. The abun- 
dance of 3 He is a crucial parameter for the intensity of 
the 2.22 MeV line emission since, though it is estimated 
to be ss 10 5 less abundant than H, it has a cross section 
for neutron capture 1.61 x 10 4 times that of H. Conse- 
quently, this isotope is able to remove neutrons from the 
atmosphere rather efficiently, thereby reducing the rate of 
their capture by protons. The density and the composition 
also determine the mean free path of 2.22 MeV photons 
in the atmosphere. 

Simulations of neutrons with an energy E n impinging 
on the atmosphere at a zenith angle ip have also been 
performed. For each assumption of E n and ip, a depth dis- 
tribution of neutron capture by H (dfd(z; E n ,ip)/dz) and 
an energy and zenithal distribution of neutrons leaving 
the atmosphere (dF(E,6; E n ,ip)/(dEd9)) are computed. 
These distributions are normalized to the number of im- 
pinging neutrons. Figure ^ shows the depth distribution 
of the sites of neutron capture by protons for neutrons im- 
pinging with an angle of 40 degrees. For neutron energies 
lower than sa 0.01 MeV, the depth distribution of neutron 
capture sites and the energy and zenithal distribution of 
neutrons leaving the atmosphere are rather constant (with 
respect to energy) for a given initial zenith angle. In fact, 
at these energies, the neutrons are quickly thcrmalized in 
the vicinity of the surface and all follow the same thermal 
diffusion. 

Depending on their energy and direction, some of the 
neutrons that are ejected out of the atmosphere get decel- 
erated gravitationally and fall back onto the secondary's 
surface; some of them decay in the process. This frac- 
tion {F r {E n ,ip)) of neutrons returning to the atmosphere 
depends on the energy and zenith angle of the imping- 
ing neutrons. Indeed, the larger the neutrons' energy, 
the deeper they can go, and the smaller the zenith an- 
gle of the impinging neutrons the deeper they will also 
go. Consequently, the probability for a neutron to leave 
the atmosphere decreases with increasing energy and de- 
creasing zenith angle. The fraction F r (E n ,ip) is calcu- 
lated using Equation (19) of Hua & Lingenfelter (1987), 
which determines the time required for a neutron to re- 
turn to the atmosphere as a function of its energy and 
zenith angle. This time is used to determine the fraction 
of neutrons decaying in flight before they return to the 
atmosphere. The spectrum of the escaping and return- 
ing neutrons has a cut-off value below the escape energy 
^ E {esc) = GMm^ w 2 fceV^g^) because even if a neutron 
has a kinetic energy a little less than the escaping energy, 
the time spent to return to the atmosphere can be larger 
than its lifetime. 

The angular distribution of the returning neutrons is 
very close to an isotropic distribution. Simulations have 
thus been performed to derive the depth distribution of 
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Fig. 2. Distributions of the neutron capture sites in the 
secondary's atmosphere for several impinging neutron en- 
ergies. 

2.22 MeV creation sites (df r (z)/dz) due to returning neu- 
trons. As previously stated, this distribution is flat for 
low neutron energies. A fraction F r p of these returning 
neutrons can again leave the secondary. However they all 
return to the atmosphere, since their energies are smaller 
than the escaping energy, and they can then again un- 
dergo one of the following: generate 2.22 MeV photons, 
decay, be captured by 3 He, or leave again the secondary 
atmosphere (with a fraction given by F r _o)- 

The total depth distribution for the radiative capture 
of neutrons is obtained by adding to the "direct" depth 
distribution (without taking into account the returning 
neutrons) the returning-neutrons component. This com- 
ponent is obtained by calculating the sum of the geomet- 
ric series of the fraction F r ^Q. The total depth distribution 
of 2.22 MeV creation site is then: 



df(z-,E n ,ip)_df d (z;E n ,ilj) df r (z) F r (E n ,ip) 



dz 



dz 



dz (1-F r , ) 



(12) 



The angular distribution of the escaping 2.22 MeV pho- 
tons is computed using the total depth distribution of 2.22 
MeV photon creation sites. Let £ and \ be the zenithal and 




Zenith angle of impinging n: 0" 270 

Fig. 3. Angular distributions of the 2.22 MeV photons 
emitted from the secondary atmosphere for several im- 
pinging neutron energies. 

azimutal direction angles of the 2.22 MeV photon. The an- 
gular distribution of escaping 2.22 MeV photons induced 
by a neutron of energy E n impinging on the atmosphere 
with a zenith angle of tp is: 



de 2 .2(€,X; E n,ip) _ J_ 
dfl 4tt 



df(z;E n ,jj) 
dz 



■ g A 2.2 c 



^ dz (13) 



with A2.2 being the mean free path of 2.22 MeV photons 
in the secondary's atmosphere. This angular distribution, 
expressed in photons steradian -1 neutron -1 , is indepen- 
dent of the azimuth. Figure || shows an example of 2.22 
MeV photons angular distribution for neutrons impinging 
with various E n 's at a zenith angle of 0°. 

The total 2.22 MeV intensity depends on the direction 
of observation with respect to the binary system frame. 
Indeed, the 2.22 MeV photons are emitted only from the 
surface of the secondary that is irradiated by neutrons. 
Therefore, the observable gamma-ray flux comes only from 
the fraction of the secondary area that is irradiated by 
neutrons and is visible for the observer (see Fig. ^). So 
we can expect to observe a periodic 2.2 MeV line flux due 
to the rotation of the binary system. The total 2.22 MeV 
intensity l2.2(&,8) (photons s -1 sr -1 ) in a direction (a, 6) 
can be obtained by estimating the integral: 



h.2{ct,5) 



Nn 
47T 



dF(E) de 2 .2^x;E^) 



D J E c 



dE 



dn 



dndE (14) 



where dF(E)/dE is the energy distribution of neutrons 
reaching the secondary, N n is the neutron rate (in n s -1 ) 
emitted by the accretion disk, <p and 9 the direction of 
the emitted neutrons (see Fig. ||) and D is the integration 
domain (<p, 8 £ D) defined as the intersection of the irra- 
diated area (delimited in a cone of angle ip m ax) and the 
visible area in the direction a and i5. In this equation, £ is 
a function of 95, 9, a and 8. The zenith angle ip depends 
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observer direction 




Fig. 4. Schematic view of the irradiation of the secondary 
star by neutrons. The binary system frame is represented 
by the dashed axes. 

Table 6. Characteristics of the secondary star for the two 
models. 





Model 1 


Model 2 


Radius (Rq) 


1 


3.6 


Mass(M ) 


1 


10 


Density (g/cm 3 ) 


ten 4 


10~ 5 


Temperature (K) 


5000 


10000 


[H] 


0.9 


0.1 


[He] 


0.1 


0.9 



only on ip. The integration domain depends not only on 
the direction of emission but also on the radius of the sec- 
ondary and on its distance with respect to the neutron 
source. 

4. Results 

We have considered two models for the composition and 
the geometrical characteristics of the secondary star (see 
Table ||) . The values shown in the table have been chosen 
so as to compare neutron capture efficiency in Helium- 
enriched and non-enriched ("normal") atmospheres. Cal- 
culations of the 2.22 MeV emissivity with specific, accu- 
rate secondary star parameters of known binary systems 
will be presented in a separate publication (Jean et al. 
2001, in preparation). 

Calculations based on the method presented in the pre- 
vious section have been performed in the aim of deriving 
the expected flux in the 2.22 MeV line from the binary sys- 
tem. The energy distribution of neutrons was assumed to 
be a Maxwellian with kT = 10 MeV, truncated at 5 MeV 
(E c in equation [l^) because the gravitationnal potential of 
the primary compact object retains the neutrons of lower 
energy (see Sect. 2). Since we are dealing with close bi- 
nary systems, the separation is less than 40 R Q , which 



is the mean free path of 5 MeV neutrons in the vacuum. 
Therefore, the energy distribution does not necessitate an 
additionnal cut to account for the loss of neutrons by de- 
cay during their trip to the secondary star. 

4-1- Fraction of escaping 2.22 MeV photons 

The fraction of escaping 2.22 MeV radiation ($2.2) is de- 
fined as the ratio of 2.22 MeV photons that reach the sec- 
ondary star surface (outward) without scattering to the 
total number of incident neutrons. This fraction is calcu- 
lated by integrating the angular distribution of escaping 
2.22 MeV photons (see equation 4) over all the possible 
angles. The 2.22 MeV photons produced by the return- 
ing neutrons are also taken into account. The fraction ^2.2 
depends on the zenith angles and energies of incident neu- 
trons. It gives an idea of the efficiency of neutrons in pro- 
ducing an observable 2.22 MeV line. Figures |^ show the 
escaping 2.22 MeV fraction for the two models of the sec- 
ondary star. The difference in the neutron capture rate 
between the two models is principally due to the differ- 
ence in the composition of the atmosphere. Indeed the 
abundance of 3 He in the second model is such that a large 
fraction of neutrons are captured by this nucleus and thus 
lost. Consequently, the neutrons capture rate decreases 
by a factor « 17 between the two models. In Model-2 the 
2.22 MeV photons are emitted deeper in the atmosphere, 
and so their shorter mean free path lowers the fraction 
of escaping photons by a factor 1.5 compared to Model- 
1. Finally, the gravitational escaping energy of neutrons 
in Model-2 is slighty lower than in Model-1, increasing 
the fraction of escaping neutrons and further reducing the 
2.22 MeV emissivity in Model-2. All these effects lead to a 
decrease of /2.2 by a factor w 30 for Model-2 with respect 
to Model-1 at low E„ . 

The fraction of escaping 2.22 MeV photons depends 
not only on the neutron capture rate but also on the 
transmission of 2.22 MeV photons in the atmosphere. In 
Model-1 (see Fig. ||a), where the neutron energy is lower 
than the escaping energy, the neutrons thermalize quickly 
and stay in the upper layer of the atmosphere. If they 
leave the atmosphere, they cannot escape the secondary 
but a fraction of them decay before returning in the atmo- 
sphere. The fraction of escaping neutrons increases with 
the zenith angle. Indeed, high-zenith-angle incident neu- 
trons can more easily escape since they scatter close to the 
surface. This is not the case for low-zenith-angle incident 
neutrons that go deeper in the atmosphere: even if they 
scatter toward the surface, their energies are rarely above 
the escape value. This explains why the /2.2 fraction de- 
creases with increasing zenith angle below 5 MeV. Above 
5 MeV, neutrons penetrate and are captured deeper in the 
atmosphere. The fraction of escaping neutrons is reduced, 
as well as the transmission of 2.22 MeV photons. However, 
for large zenith angles, the neutron captures happen closer 
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Fig. 5. Fraction of escaping 2.2 MeV photons as a function 
of the impinging neutron energies for several impinging 
zenith angle. The two models are presented separately. 



to the surface in zones allowing the 2.22 MeV photons to 
escape more easily from the star. 

The variation of the 2.22 MeV photon escaping frac- 
tion as a function of the angle and the energy of impinging 



neutrons in Model-2 is similar to that in Model- 1. How- 
ever, a resonance in the elastic scattering of neutrons with 
4 He at 1.1 MeV is at the origin of the drop around 1 MeV 
(see Fig. |b). At this resonance, the angular distribution 
of 4 He-scattered neutrons is peaked for backward direc- 
tions Q Consequently, around 1 MeV, incident neutrons 
are more likely to backscatter and to be ejected out of the 
atmosphere with kinetic energies too high to allow them 
to return. This effect reduces the number of MeV neutrons 
that can thermalize and produce 2.2 MeV photons, hence 
the trough seen in Fig. ||b. 

4-2. Flux versus model and irradiation geometry 

The distance between the accretion disk (the source of 
neutrons) and the center of the secondary star, and the 
angles a & S, which specify the direction of Earth (see 
Fig. H), are the basic parameters of this calculation. Since 
the binary system is rotating, a varies with time and the 
observed flux becomes periodic. Intensities have been cal- 
culated as a function of the direction of emission for dif- 
ferent values of the distance between the primary and sec- 
ondary stars. Assuming parameters of the binary system 
(separation and masses), it has been possible to derive 
light curves of the 2.22 MeV line emission, depending on 
5. However, the range of interesting separations is limited 
since the secondary radius needs to fill its Roche lobe in 
order to allow the accretion of its matter by the compact 
object. Using a study by Paczynski (1971) that calculated 
the critical radius of a star in a binary system that allows 
the outward flow of its matter, we estimated the upper 
limit of the separation to be between 3.5-6 R Q for a 
secondary star of 1 M Q (Model-1) and 7.5 - 11 R Q for 
a secondary star of 10 M Q (Model-2), providing compact 
object masses ranging from 3 to 20 M . With this type of 
close binary systems, the period of rotation is less than a 
day. 

Figure H shows phasograms of the 2.22 MeV intensity 
as a function of the separation for Model-1. As expected, 
the intensity is maximum when the neutron-irradiated 
area of the secondary star atmosphere faces the observer. 
The maximum value of the intensity decreases with in- 
creasing separation values because the neutrons flux im- 
pinging on the atmosphere decreases with increasing dis- 
tance between the accretion disk and the secondary. In 
Fig. ^, the intensity is equal to zero at a phase value of 
0.5 only for a separation of 1.5 R , because for this value 
of S the irradiated area is hidden by the star. For larger 
separations, the irradiated area becomes larger and a frac- 
tion of this area is always visible by the observer. 

The mean and the root-mean-square (RMS) fluxes are 
commonly used for the analysis of periodic emissions. 
They have been estimated for our 2.22 MeV intensity as 
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Fig. 6. Example of phasogram of the 2.2 MeV emission for 
the Model-1 with a compact object mass of 3 Mq. S is the 
angle between the binary system plan and the direction of 
the observer. D is the distance between the accretion disk 
and the secondary star in solar radius unit. The period of 
the binary system change from 0.1 to 0.4 days. 

a function of the separation and the direction of the ob- 
server. Both the mean and the RMS intensities decrease 
with increasing separation values as a power law with slope 
around -2 (between -1.8 and -2.3). Figures fj] and |§] show 
the variation of the mean and RMS intensities as a func- 
tion of the binary separation for the two models. 

The mean intensity does not vary significantly with 
the angle of the Observer's direction with respect to the 
binary system plane. The RMS intensity decreases, how- 
ever, showing a reduction of the modulation of the visible 
photon flux with increasing d. The modulation obviously 
disappears for a direction of observation of the binary sys- 
tem perpendicular to the plane of the latter (<5=90°) . 

The flux at 2.22 MeV can be derived using the formula: 

ci in— 5 ^™ -*2.2 ,_2 -1 -2 /-i r\ 

-F2.2 = 10 — -375 t — -3 5 r d k s cm (15) 

10 4U ns^ 1 10~ 2 n^ 1 sr^ 1 p 

where R n is the rate of neutrons emitted by the accretion 
disk, I2.2 is the normalized intensity and dk pc the distance 
of the binary system to the observer in kpc. For Model-1, 
the flux at 1 kpc is found to be in the range 10 -7 - 10 -5 
photons s _1 cm~ 2 depending on the accretion disk model. 
For Model-2, it is found to be in the range 10 -8 to 10~ 6 
photons s" 1 cm~ 2 depending on the accretion disk model. 
Combining the results of R n from Sect. 2 (Tables la,b to 
5a, b) with the calculations in this section, we obtain the 
final fluxes F 2 .2, which we present succinctly in Tables 7 
and 8. The separations have been chosen to be 2 R and 8 
Rq for Model-1 and Model-2 respectively, since for these 
values the mean fluxes do not change significantly with 
the binary system inclination (see Figs. and |J). 

The rotation of the binary system leads to a shift in 
the ccntroid of the 2.22 MeV line. This Doppler shift 
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Fig. 7. Normalized mean (a) and RMS (b) intensities for 
Model-1 as a function of the separation of the binary. Sev- 
eral values of 5 are shown. 



changes with the phase and depends on the observer's 
direction (6). Using classical relations in binary systems 
(e.g. Franck, King & Raine, 1992) we derive a relation 
(Equation [If]) for the spectral shift of the 2.22 MeV line 
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which depends on the masses of the compact object and 
the secondary (M c and M s , respectively, in solar masses), 
the separation between the two objects (d s in solar radii), 



Table 7. Mean 2.22 MeV flux vs. Accretion disk models 
(at 1 kpc for a separation of 2 R and an initial 90% H 
and 10% He Composition) 



Accretion model M (Mq jyr) 


a = 0.3 


a = 0.1 




10 -io 


1.1 x 10~ 8 


5.5 x 10" 8 


ADAF 


10~ 9 


4.3 x 10~ 7 


1.4 x 10~ 6 




10~ 8 


9.1 x 10 


1.5 x 10~ 5 




10 -io 


1.6 x 10" 8 


1.3 x 10~ 7 


ADIOS 


10 


1.5 x 10" 6 


7.8 x 10 -6 




10~ 8 


8.3 x 10~ 5 


1.1 x 10 -4 




10 -io 


3.5 x 10 -22 


4.8 x 10 -12 


SLE 


10" 9 


4.1 x 10" 9 


2.8 x 10~ 7 




10 -8 


4.8 x 10" 6 


1.8 x 10~ 5 






kT e = 0.5MeV 


kT e = O.lMeV 




10 -io 


3.1 x 10" 7 


4.7 x 10" 8 


UIT - 30MeV 


10" 9 


2.3 x 10" 6 


4.2 x 10~ 7 




10" 8 


1.6 x 10" 5 


4.9 x 10~ 6 




10 -io 


1.9 x 10~ 7 


2.0 x 10~ 8 


UIT - lOMeV 


10" 9 


1.6 x 10" 6 


2.0 x 10~ 7 




10" 8 


1.1 x 10" 5 


1.7 x 10~ 6 



the phase (4>) , and the angle between the observer's direc- 
tion and the binary system (5). 



AE = 3.2 



y/d s (M c + M s ) 



cosS sirup keV. 



(16) 



The estimated spectral shift is mesurable with SPI, the 
spectrometer onboard the INTEGRAL satellite. Using 
Equation [l6] and the intensity variation as a function of 
the phase (as presented in Figure ||) we have obtained the 
2.22 MeV line profiles for Model-1 (Figures ^) with a com- 
pact companion of 3 Mq and for two values of the binary 
system inclination^]. The line is broadened and double- 
peaked due to the rotation of the secondary star. Figures 
|To| show the shape of the line as it could be measured 
by SPI (assuming a spectral resolution of 2.9 keV at 2.22 
MeV) . For an inclination of 80° the double-peak shape is 
clearly visible, whereas for a 10° inclination that is not 
the case. 



5. Conclusions 

We have presented a theoretical investigation of the 2.22 
MeV line emitted from X-ray binary systems neutrons are 
captured in the secondary's atmosphere. These neutrons 
are produced in the accretion disk by nuclear spallation of 



2 The inclination i is the angle between the binary system 
rotation axis with respect to the observer direction, i.e. i — 
90° - S. 
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Table 8. Mean 2.22 MeV flux vs. Accretion disk models 
(at 1 kpc for a separation of 8 R0 and an initial 10% H 
and 90% He Composition) 



Accretion model M (Mq /yr) 


a = 0.3 


a = 0.1 




10- 10 


1.0 x 10 -9 


4.0 x 10" 9 


ADAF 


10 9 


3.9 x 10~ 8 


9.2 x 10~ 8 




10~ 8 


6.0 x 10~ 7 


1.1 x 10 -6 




10- 10 


2.1 x 10~ 9 


1.2 x 10~ 8 


ADIOS 


10 


1.3 x 10 -7 


5.2 x 10~ 7 




10- 8 


5.4 x 10~ 6 


1.5 x 10~ 5 




10 -io 


9.7 x 10" 23 


1.3 x 10~ 12 


SLE 


1Q -9 


7.5 x 10" 10 


2.7 x 10~ 8 




10" 8 


4.9 x 10 -7 


2.0 x 10" 6 






kT e = 0.5MeV 


kT e = O.lMeV 




10 -io 


1.1 x 10" 8 


1.5 x 10" 9 


UIT - 30MeV 


10 -9 


9.0 x 10" 8 


1.5 x 10~ 8 




10~ 8 


6.6 x 10~ 7 


4.6 x 10~ 7 




1Q -io 


6.0 x 10" 9 


6.1 x 10~ 10 


UIT - lOMeV 


10~ 9 


5.4 x 10" 8 


7.2 x 10~ 9 




10~ 8 


4.1 x 10~ 7 


1.2 x 10~ 7 



Helium nuclei. The rate of production of neutrons depends 
on various parameters set by the accretion disk model, 
i.e. the temperature, density, and composition of the ac- 
creted material. Calculations of the line intensity normal- 
ized to the neutron production rate have been performed 
for two simple secondary models and for several X-ray bi- 
nary geometries. 

We showed that, due to the rotation of the secondary, 
the intensity of the 2.22 MeV radiation is periodic, and the 
line centroid is shifted by the Doppler effect. The mean 
intensity does not vary significantly with the direction of 
observation with respect to the binary system frame. 

According to our model, the spectral (Figures || and 
|l0| ) and temporal (Figures |[ and ||) analyses of the 
2.22 MeV line flux would provide valuable insights into 
the characteristics of the binary system (separation, in- 
clination with respect to the observer, masses, accretion 
rate and disk structure, etc.). If one further knows the 
X-ray luminosity, the geometry of the X-ray binary sys- 
tem, and its distance from Earth, then a measure of a 2.22 
MeV line flux from it can set constraints on the neutron 
production rate and consequently on the accretion disk 
models (ADAF, ADIOS, SLE, etc.). 

At 2.22 MeV the narrow 7-ray line sensitivity of the 
SPI spectrometer of the INTEGRAL mission is expected 
to be 7-10 xlO -6 photons s" 1 cm for an observa- 
tion time of 10 6 seconds (Jean et al, 1999). The 2.22 
MeV fluxes estimated for the simple secondary-star mod- 
els and for the various accretion disk models used here are 
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Fig. 9. Shape of the neutron capture line for Model-1. 
The inclination of the binary system is 10° (a) and 80° 
(b). Several values of the binary separation are shown. 

found to be measurable by the spectrometer SPI in some 
cases: mostly if the accretion rate is large enough (M « 
10 _8 M©/yr) and the viscosity is significant (a £ 0.2). In 
case of a detection, the SPI spectral resolution (ss 3 keV 
at 2.22 MeV) would allow the measurement of the broad- 
ening of the line, which is due to the rotation of the sec- 
ondary, if the separation and the inclination of the X-ray 
binary are not too large. 

In a future work, the 2.22 MeV line flux and shape 
will be estimated for known X-ray binary systems (e.g. 
A0620-00, etc.) using all accurate, available information 
on their characteristics (distance, separation, composition, 
etc.). The calculated fluxes will be compared with upper- 
limit fluxes obtained with COMPTEL (Van Dijk, 1996), 
and correlation with other gamma-ray lines will be in- 
vestigated. Indeed, the neutrons that irradiate the sec- 
ondary star can also produce Be and Li isotopes in its 
atmosphere. This process was proposed by Guessoum & 
Kazanas (1999) to explain the overabundance of Li in 
some soft X-ray transients. Both 7 Li and 7 Be emit gamma- 
ray lines (at 0.478 and 0.429 MeV), and 7 Be decays into 
7 Li* (with a half- life of 56 days), the latter then produces 
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Fig. 10. Estimation of the neutron capture line shape for 
Model- 1 as it could be measured by the spectrometer SPI 
of INTEGRAL. The inclination of the binary system is 10° 
(a) and 80° (b). Several values of the binary separation are 
shown. 



0.478 MeV radiation upon de-exciting; such a delayed 
gamma-ray emission, although weak, would be a clear sig- 
nature of the process. We can thus expect a certain cor- 
relation of the 2.22 MeV emission with these gamma-ray 
lines, providing supplementary insights into the physics of 
the binary system. 
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